# 导入sympy
from sympy import *
# 初始化打印格式，让矩阵以常见方式输出
init_printing(use_unicode=True)
# 注意矩阵的输入格式
A = Matrix([[1,5,-1,-1],[1,-2,1,-3],[3,8,-1,1],[1,-9,3,7]])
# 输出矩阵
print(A)
# 求矩阵A的行数和列数
print("矩阵A的行数和列数为：",A.shape)
# 求矩阵A的转置矩阵
print("矩阵A的转置矩阵为：\n",A.T)
# 求矩阵A的3次幂
print("矩阵A的3次幂为：\n",A**3)
# 求矩阵A的行最简阶梯型矩阵
print("矩阵A的行最简阶梯型矩阵：\n",A.rref())
# 求矩阵A的秩
print("矩阵A的秩为：",A.rank())
# 求矩阵A的逆矩阵
print("矩阵A的逆矩阵为：\n",A.inv())
print("矩阵A的逆矩阵为：\n",A**(-1))
# 报错说明没有逆矩阵